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INTRODUCTION 

Massive fermions have long been a problem in gauge theories. Unification of elec- 
tromagnetic and weak forces was once hindered by the fact that the introduction of 
mass terms broke the gauge invariance of the theory. This problem was solved by the 
introduction of the Higgs field. Spontaneous breakdown of the SU(2) x U(l) symme- 
try then takes place. The gauge bosons gain mass and the masses for the fermions are 
generated through their Yukawa interaction with this Higgs field. However, there has 
been a widespread dissatisfaction with this mechanism since the masses are not pre- 
dictable. Rather, they must be fixed by experiment. Studying the non-perturbative 
behaviour of gauge theories provides an alternative. If the interactions are strong 
enough, they are capable of generating masses for the particles dynamically even if 
they start with zero bare mass. Moreover, experiment tells us that the top quark is 
very heavy and so the Yukawa coupling gt for top-Higgs interaction is 0(1). Then 
one naturally expects that non-perturbative effects become important. Indeed, it has 
been suggested [1] that the top quark may acquire mass non-perturbatively through 
four-fermion interactions, and the Higgs can then be viewed as the condensate of the 
top and the antitop. However, in an attempt to include the effects of gauge boson ex- 
change term, one loses gauge invariance of the physical quantities. Of course, physical 
quantities must be gauge independent. This motivates the study of how to achieve 
this in non-perturbative calculations. Quenched QED provides a toy model in which 
to study this problem, as we discuss. 



DYSON-SCHWINGER EQUATIONS 

Our starting point is the set of Dyson- Schwinger equations. These are an infinite 
system of coupled equations for all the Green's functions, which are non-perturbative 
in nature. Their structure is such that the 1-point function is related to the 2-point 
function, the 2-point function is related to the 3-point function, etc. ad infinitum. As 
it is impossible to solve the complete set of equations, one has to truncate this infinite 
tower in a physically acceptable way to reduce them to something that is soluble. A 
familiar way to do this is perturbation theory. However, if one wishes to generate 
masses for particles, a non-perturbative way has to be sought. 

To see how to do this consider two of the Dyson- Schwinger equations, one for the 
fermion propagator, and the other for the photon propagator. These are shown below 
diagrammatically together with their corresponding mathematical expressions: 
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FIG. 1. Dyson-Schwinger equation for fermion propagator. 
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FIG. 2. Dyson-Schwinger equation for photon propagator. 
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where the quantities with the superscript '0' are bare quantities, and the others are 
full ones. Quenched QED corresponds to making the assumption that the full photon 
propagator can be replaced by its bare counterpart. This limit is achieved by regarding 
Nf as a mathematical parameter, which is set equal to zero. As an example, to begin 
with, we make a further simplification by replacing the full vertex by the bare one. 
Eq. (1) then reduces to: 
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FIG. 3. Rainbow approximation. 
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in what is known as the rainbow approximation, where 



Spik) 

sUk) 
KM) 



|(f-M(A;2) 
1 



— mo 

1 

~2 



Q 



9.^ + {^-l)'-f 



Eq. (3) is a matrix equation which corresponds to two equations in A4 and F. We 
can project out equations for these by taking the trace of Eq. (3) having multiplied 
by 1^ and 1 in turn to obtain: 
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where as usual a — e^/Ai:. On carrying out the angular integrations, and putting the 



bare mass equal to zero, we have 
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where A is the ultra-violet momentum cutoff. It is easiest to solve these equations 
in the Landau gauge where they decouple. F{p'^) is obviously 1. Moreover, there 
is a non-trivial solution [2] for the mass function Ai for the coupling larger than a 
critical value of etc = 7r/3. This is best illustrated by plotting the Euclidean mass 
M — M.{M'^) as a function of a, as found by Curtis and Pennington [3]: 




FIG. 4. Euclidean mass, M = Ai{M'^) dynamically generated in the rainbow approximation 
as a function of the coupling a in three different gauges: Landau = 0) •, Fcynman = 1) 
A, and Yennie = 3) o gauges. 

Note that = is always a solution to Eq. (5). However, beyond the critical 
value of the coupling, the non-zero solution bifurcates away from the trivial solution. 
This is in complete contrast with the perturbation theory, where, even if we perform 
an all orders resummation using the Renormalization Group Equation, we end up 



with a result of the following form, 
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and the field remains massless to all orders if we start with a zero bare mass, mo = 0. 

In contrast, non-perturbative dynamics is able to generate masses for particles 
even if they have zero bare mass. However, there are problems. As the critical 
coupling corresponds to a change of phase, we expect it to be independent of the 
gauge parameter. But when one solves the Eqs. (4) and (5) for different gauges, one 
finds that this is not the depicted in Fig. 4. However, it is not difficult to trace 

the root of this problem. The full vertex of Eq. (1) has to satisfy the Ward-Takahashi 
identity for the fermion propagator to ensure its gauge covariance. However, the bare 
vertex that was used in Eq. (3) does not obey this identity. Therefore, one should 
not expect physical outputs to be gauge independent when the input is not. 

THE VERTEX 

We expect that any reasonable ansatz for the vertex should fulfill the following 
requirements: 

• It must satisfy the Ward-Takahashi Identity in all gauges. 



• It must ensure that the fermion propagator of Eq. (1) is multiplicatively 
renormalizable . 

• It must result in a critical coupling, at which mass is generated dynam- 
ically, that is gauge independent. 

• It must be free of any kinematic singularities, i.e. it should have a 
unique limit when k'^ — > . 

• It must have the same transformation properties as the bare vertex 7^ 
under C and P. 
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Keeping in mind the form of the Ward-Takahashi identity, one can spht the full vertex 
into two components, longitudinal and transverse: 

r^{k,p) = r^{k,p) + ri^{k,p) , (6) 

where, the transverse part of the vertex is defined by: 

q^r!^{k,p) = 0. (7) 

The Ward-Takahashi identity uniquely fixes the longitudinal part of the vertex, as 
shown by Ball and Chiu [4], to be 
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However, the transverse part remains arbitrary. Ball and Chiu [4] enumerated a basis 
of eight independent tensors in terms of which the most general form for the transverse 
part of the vertex can be written: 

mk,p)^j2n{e,p',q')Tt{k,p) . (9) 
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We list here only those four tensors which we shall need later: 

T,^{k,p) = {p^{k.q)-k^{p.q)){f^+^) 

T^Kp) = q^l'-q'i 

mk,p) = rik'-p')-ik+pm-l() 

T^{k,p) = -Yp'k^'a.p+jf ^-k^ ^ , (10) 

with (7fj,i, — The simplest choice is to take the transverse part to be zero. 

But Curtis and Pennington [5] showed that if we take the transverse part of the vertex 



to be zero, the fermion propagator is no longer multiplicatively renormalizable. They 
suggested the following transverse part of the vertex satisfying this requirement. 
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where, d{k'^,p'^) — k^ for k"^ ^ p^. d{k'^,p'^) must be symmetric in k and p and free of 
kinematic singularities leading to the proposal: 
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The vertex specified by Eqs. (8-12) will be referred to as the CP-vertex [5]. Curtis 
and Pennington solved the coupled equations for F and A4 from Eq. (1), using 
this ansatz. They found that the gauge-dependence of the critical coupling at which 
the non-perturbative behaviour bifurcates away from the perturbative one reduces 
considerably, as seen by comparing Figs. 4 and 5. 
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FIG. 5. Euclidean mass, M = AI(M^) dynamically generated using the CP-vertex as a 
function of the coupling a in three different gauges: Landau (^ = 0) •, Feynman (^ = 1) A, 
and Ycnnie (^^ = 3) o gauges. This plot is to be compared with the rainbow approximation 
results of Fig. 4. 



BIFURCATION ANALYSIS 



To see this, Atkinson et al. [6] recently suggested a bifurcation analysis to study 
the phase change near the critical coupling. This is a precise way to locate the 
critical coupling as compared to the previous methods which rely on numerical cal- 
culations. This method amounts, in practice, simply to throwing away all terms that 
are quadratic or higher in the mass-function A4. Employing this procedure, and using 
the fact that at the critical coupling, M{p'^) ~ and -F(p^) ~ {p'^Y in Eq. (1), 

one arrives at the following equation in an arbitrary gauge: 
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There are two roots of this latter equation for s between and 1. Bifurcation occurs 
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FIG. 6. Critical coupling, Uc, as a function of the gauge parameter, ^ (solid line). The 
corresponding values for the rainbow approximation have also been shown o. 



The bifurcation point defines the critical couphng, etc- Numerically, etc = 0.933667 
in the Landau gauge. For each value of the gauge parameter, these equations can be 
solved for v, Sc and etc- The solution found by Atkinson et al [6] is displayed in Fig. 
6. For comparison, the points for the bare vertex have also been shown. One can 
see that the gauge dependence has considerably been reduced, as was seen earlier. 
However weak this variation, any gauge dependence shows that the CP vertex cannot 
be the exact choice. 

CONSTRAINTS OF MULTIPLICATIVE RENORMALIZ ABILITY 

To find a vertex that ensures the gauge independence of the critical coupling, we 
start off by making three assumptions. Firstly, we demand that a chirally-symmetric 
solution should be possible when the bare mass is zero, just as in perturbation the- 
ory. This is most easily accomplished if the sum in Eq. (9) involves just i = 2, 3, 6 
and 8. The second assumption is that the functions, Tj, multiplying the transverse 
vectors, Eq. (9), only depend on k'^ and p^, but not q'^. The third assumption is that 
the transverse part of the vertex vanishes in the Landau gauge. The motivation for 
this comes from the lowest order perturbative calculation for the transverse vertex, 
satisfied by Eq. (11). These conditions fix the of Eq. (9). Multiplicative renor- 
malizability of the wavefunction renormalization F(p^) enables us to write Tq and r 
in terms of one function Wi{x) [7] : 
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where r is the combination of Tj given by: 
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The condition of multiplicative renormalizability, i.e, -F(p^) ~ {p^Y, constrains the 
otherwise arbitrary function Wi as follows: 



It should be noted that, with the simplest choice Wi — 0, the massless CP- vertex, 
Eqs. (11,12) emerges. 

CONSTRAINTS OF GAUGE INVARIANCE 

At the bifurcation point, as stated before, multiplicative renormalizability forces 
a simple power behaviour for the mass function as well as for the wavefunction renor- 
malization. Such a multiplicatively renormalizable mass function must exist in all 
gauges. Consequently, the exponent, Sc, must be gauge independent. Moreover, dy- 
namical mass generation marks a physical phase change and so the critical coupling, 
etc, must also be gauge independent. Thus the critical values, etc, Sc, found in the 
Landau gauge must hold in all gauges. Using this physically motivated argument, the 
equation for the mass function gives T2, T3 and ts in terms of a function W2{x) [8]: 
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and the function W2 is constrained, by the gauge invariance of the mass function and 
the critical couphng, to obey the following integral equation, 

dx 
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at the critical coupling a — ac- In order to make sure that none of the functions Tj 
has kinematic singularities as /c^ P^ ■, M^i and 1^2 should also satisfy the following 
conditions: 

W^2(1) + 2W^^(1) = 2i{v-s-\) . 



This defines the construction of the full vertex via Eqs. (6-10) [8]. 



CONCLUSIONS 



Above we have presented a truncation of the fermion Schwinger-Dyson equa- 
tion for the quenched QED, which respects the key properties of the theory. We 
have constructed a non-perturbative vertex in terms of the constrained functions 
Wi{x){i = 1,2). It satisfies the Ward-Takahashi identity, ensures the fermion prop- 
agator is multiphcatively renormahzable, agrees with one loop perturbation theory 
for large momenta and enforces a gauge independent chiral symmetry breaking phase 
transition. This study motivates the need for a realistic investigation of tt conden- 
sates as the source of the electroweak symmetry breaking. Including the four fermion 
interaction, the Dyson- Schwinger equation for the fermion propagator becomes: 



FIG. 7. Dyson- Schwinger equation for the fermion propagator, including the four fermion 
interaction term. 

We need to solve this equation in a gauge invariant way. The study of quenched 
QED presented here suggests that a proper choice of the vertex can guarantee the 
gauge independence of the physical observables. However, a reahstic calculation, of 
course, requires the unquenching of the theory which complicates the problem signif- 
icantly. The fermion-boson vertex (in particular its transverse part) will intimately 
depend on the photon renormalization function in a non-perturbative way not yet 
understood. The discussion for quenched QED presented here provides the starting 
point for such an investigation of full QED. 
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